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Aims. We measure the mass and kinematics of the massive molecular torus in the planetary nebula NGC 6302. The nebula is the proto-typical 
butterfly nebula. The origin of the wing-like morphology is disputed: determining the mass-loss history of the confining torus is an important 
step in understanding the formation of this structure. 

Methods. We performed submillimeter observations with JCMT and the SMA interferometer The continuum emission as well as the J=2-l 
and 3-2 transitions of '-CO and '^CO are analysed at arcsecond resolution. 

Results. The CO emission indicates amass of the torus of ~ 2 Mo + IMq. The '^CO and ''CO emission matches the dark lane seen in absorption 
in the Ho- image of the object. The CO torus is expanding with a velocity of ~ 8 km s"', centred at Vi,, = -31.5 km s"'. The size and expansion 
velocity of the torus indicates that the torus was ejected from ~ 7500 yr to 2900 yr ago, with a mass-loss rate of 5 x lO"'' Mq yr"' . We also see 
a ballistic component in the CO images with a velocity gradient of 140 km s"' pc"' . 

Conclusions. The derived mass-loss history of the torus favours binary interaction as the cause of the ejection of the torus. We predict the 
existence of a companion with an orbital period P ^ I month. 

Key words. Stars: AGB and post-AGB; planetary nebulae: general; planetary nebulae: individual: NGC6302 



1. Introduction 

Asymptotic Giant Branch (AGB) stars lose between 50% and 
80% of their stellar mass at mass-loss rates up to 10""* Moyr"'. 
About 80% of planetary nebulae (NP), i.e. the visible part of 
the ionized gas ejected by AGB stars, show elliptical and/or 
bipolar morphologies (e.g. Manchado 1997). Bipolar FN are 
believed to form when a fast wind from the now post-AGB star 
ploughs into the earlier slow AGB wind, amplifying an initial 
equatorial density enhancement (Kwok et al. 1978, Balick et al. 
1987; Frank & Mellema 1994, Zijlstra et al. 2001). Some 15% 
of PNe show strong morphological evidence for such an equa- 
torial enhancement of cold gas around the central star (CoiTadi 
& Schwarz 1995). Detailed studies of these equatorial struc- 
tures can then give direct insights on the formation of bipo- 
lar PNe. Related to bipolar nebulae, butterfly PNe do not show 
well collimated outflows, but have thicker expanding equato- 
rial torii. The complex morphology of their outflows has been 
attributed to an interaction between the fast wind and a warped 
disk (Icke 2003). 

NGC 6302 is the proto-typical 'butterfly' nebula. It is lo- 
cated at a distance d ~ I kpc (Meaburn et al., 2005). A broad 
central absorption lane at its centre is suggestive of a massive 
circumstellar torus. The nebula is driven by one of the highest 
temperature PN central stars known (Casassus et al. 2000), al- 
though due to the high extinction of the absorption lane, it has 



not been directly detected. Previous CO observations suggest a 
mass for the torus of ~ 0.1 Mq (Gomez et al. 1989; Huggins & 
Healy 1989). However, recently Matsuura et al. (2005a) mea- 
sured a total gas mass of ~ 3 Mq from observations of the 
submillimetre dust continuum emission. This extremely high 
mass, especially when compared to the ~ 0.7 Mq estimated 
for the central star (Casassus et al. 2000), prompted us to re- 
examine the CO emission from this source and study of the 
kinematics of the circumstellar gas. 



2. Observations 

2.1. JCMT observations 

NGC 6302 was observed with the James Clark Maxwell 
Telescope in Hawaii as observing program S04AU08 

on 26 May, 18 and 19 June 2004 using the facility 230 GHz 
and 345 GHz receivers. We observed '^CO J= 2 ^ 1 and 
J= 3 — > 2, as well as the corresponding transitions of '^'CO. 
The data were taken with a velocity resolution of ~ 0.1 kms"' 



' The James Clerk Maxwell Telescope is operated by The 
Joint Astronomy Centre on behalf of the Particle Physics and 
Astronomy Research Council of the United Kingdom, the Netherlands 
Organisation for Scientific Research, and the National Research 
Council of Canada. 
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Table 1. Observed transitions 



Transition 




Frequency 


Telescope 


Resolution 


Beam 






(GHz) 




(kms-i) 


(arcsec) 


'-CO J= 2 


1 


220.398 


JCMT 


0.4 


20 








SMA 


0.3 


6.2x3.2 


'^CO J= 2 ^ 


1 


230.538 


JCMT 


0.4 


20 








SMA 


4.4 


6.2 X 3.2 


'^CO J= 3 ^ 


2 


345.796 


JCMT 


0.5 


13 


"CO J= 3 ^ 


2 


330.587 


JCMT 


0.5 


13 



for all the transitions, but was rebinned to 0.4 km s ' for anal- 
ysis. 

We raster-mapped the region around NGC 6302 in the 
J= 2 ^ 1 and J= 3 ^ 2 transitions of '^CO. Extended emis- 
sion components were seen at some velocities, but the main 
emission component corresponding to NGC 6302 was found 
to be unresolved at 230 GHz and only slightly resolved at 
345 GHz. We subsequently took '^CO spectra pointing to- 
wards the '^CO J= 3 ^ 2 emission peak position, taken to be 
a = 17hl3"43.9^ 5 = -37°06'11" (J2000). 

All observations were made by position-switching to an off 
position (0",-h300") from this central position. This off posi- 
tion was carefully checked from emission which may have con- 
taminated the source observations, but no emission with a peak 
temperature of greater than 0.2 K (with a line width of 1 km/s 
or more) was found at the frequency of either the J= 3 — > 2 
or J= 2 — > 1 transitions. The telescope pointing was regularly 
checked during the observations and was found to be constant 
to within +2". A linear baseline has been removed from each of 
the spectra. The main beam efficiency of the JCMT is 0.69 and 
0.63 at the frequency of the J= 2 ^ 1 and J= 3 — » 2 respec- 
tively. The observational parameters are summarised in Table 



m 



2.2. SMA observations 

NGC 6302 was observed at a wavelength of 1.3 mm with 
the SMaH on 26 June 2005. We used all 8 antennas in the 
compact configuration, with a maximum of baseline of 70 m. 
The six shortest baselines range in length from 16.5 to 32 m. 
The phase center was RA(2000)=17h''13'"44.2' Dec(2000)=- 
37°06'15.9". The central frequency of the upper side band was 
set to 230.538 GHz with a total bandwidth of ~ 2 GHz (i.e. cov- 
ering the frequency range ~ 229.5 to ~ 231.5 GHz). The lower 
side band was centered at ~ 220.5 GHz (i.e. covering ~ 219.5 
to ~ 221.5 GHz). The correlator was configured to provide a 
velocity resolution of ~ 0.26 kms ' for the CO J= 2 — > 1 line 
and ~ 4.4 kms ' for the '^CO J= 2 — » 1 line. Weather condi- 
tions were good for 230 GHz observations with a relative hu- 
midity of 20%, atmospheric opacity at zenith T225 - 0.2 and 

Tsys.DSB ~ 250 K. 




^ The Submillimeter AiTay is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics, and is funded by the 
Smithsonian Institution and the Academia Sinica. 



Fig. 1. Ho- image from Matsuura et al. (2005a) (colour scale) 
overlaid with our SMA 1.3 mm continuum map (white con- 
tours). The contours are 5% and then from 10 to 90% in steps 
of 20% of the 1.3 mm continuum emission peak. The emission 
peak is 0.62 Jy/beam. 

The phase and amplitude were calibrated by observing 
SgrA* for 5 minutes around every 20 minutes of integration 
on NGC 6302. The angular separation between SgrA* and 
NGC 6302 is ~ 11°. The passband calibration was done us- 
ing Callisto, 3C279 and 3C454.3. The flux calibration was per- 
formed using Callisto assuming a brightness temperature of 
120 K at 229.45 GHz. The absolute astrometric positional ac- 
curacy in the observations is found to be better than 0.3"based 
on observations of the quasars 191 1-201 and 1924-292. 

The visibility data were reduced using the SMA version 
of the Miriad software package. Applying natural weighting, 
the final maps have a synthesized beam of 6.2"x3.2"with 
P.A.=6.9°. The rms noise in the continuum image obtained by 
averaging the LSB and USB data is cr = 5.9 mJy. The rms noise 
in the CO channel maps, at 1 kms ' resolution, is ~ 125 mJy. 

3. Continuum emission 

The SMA 1.3 mm continuum map is shown in Fig.[T] To avoid 
molecular line contamination of our continuum data, the CO 
emission, on which the the correlator bandwidth was centred, 
has been removed. This continuum map reveals a single source 
elongated in the North-South direction. The angular resolution 
provided by the SMA spatially resolves the structure. The po- 
sition, size and flux are given in Table |2l together with the inte- 
grated flux measured inside the 5cr contour (~ 30 mJy/beam). 

The 1.3 mm continuum flux contains contributions from 
dust emission and thermal free-free emission. The free-free 
contribution is ~ 50% of the total flux (see Fig. 13 of Matsuura 
et al. 2005a), leaving 0.65 Jy from the dust continuum emis- 
sion. Assuming a mean dust temperature ranging from 30 K 
to 120 K (Kemper et al. 2002), a dust to H2 mass ratio of 
1 %, a typical dust opacity for protostellar disks of 2 cm^ g~' 
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(Beckwith et al. 1990), and a distance d ~ I kpc (Meaburn et 
al., 2005), we infer a total (dust plus gas) mass ranging from 
0.8 Mq to 3.9 Mq. However, the temperature range given by 
Kemper et al. (2002) corresponds to two different components, 
a cold one with = 30 to 60 K, and a hot one with = 
100 to 120 K. As discussed later (see Sect. 5.3), the bulk of the 
mass comes from the cold component, meaning that 0.8 Mq is 
rather a lower limit to the mass. Taking 60 K as the average 
temperature, we find a mass of 1.8 Mq. Uncertainties due to 
temperature and dust opacity are at least a factor of 2 in each 
way. Comparing this mass estimate to the one from Matsuura 
et al. (2005a) obtained from single dish observations, we con- 
clude that there is almost no missing dust continuum flux in our 
SMA data, as confirmed later in Sec. 5.1 through our molecu- 
lar line data set. Because the SMA is not sensitive to extended 
emission larger than approximately 20" the lack of missing flux 
confirms then the compact nature of the NGC 6302 torus. 

The peak of the millimetre map is slightly offset, i.e. ~ 2", 
from the centre of the dark lane. Two effects explain this off- 
set. First, the contribution from the free-free emission seen at 
6 cm (Matsuura et al. 2005a) shifts the 1 .3 mm peak toward the 
location of the free-free emission peak. Second, the dust ab- 
sorption is tracing only the foreground part of the dusty torus, 
whilst the 1.3 mm emission traces it entirely. A slight tilt of the 
torus results in the 1.3 mm emission appearing offset from the 
absorption. 

4. Line data 

4.1. Averaged spectrum 

The JCMT data show a compact CO structure approximately 
coincident with the continuum emission, which traces the 
molecular content. The CO spectra for all four transitions, at 
the central position, are shown in Fig.|2] The spectra show sev- 
eral broad components, with emission ranging between -15 
and -65kms using local standard of rest (LSR) velocities. 
In addition, several narrow components are seen in '^CO only, 
and are better seen in J= 2 ^ 1 than in J= 3 — > 2. These com- 
ponents are at -38, -30 and -10km s"'. The relative strengths 
of these features suggest they trace gas which has a high 
'^CO/'^CO ratio, and may be spatially extended. 

The SMA has both higher resolution, and as an interferom- 
eter it is not sensitive to very extended emission. This allows 
us to identify interstellar emission components, which will be 
absent from the SMA maps. The averaged '^CO J= 2 — > 1 ob- 
tained with the SMA is shown in Fig. [3]El Two of the nar- 
row components are absent, but the one at -38kms"' is seen, 
suggesting this gas is associated with the molecular torus. Its 
weakness at J= 3 — > 2 in the JCMT data can be attributed to an 
off-centre position combined with the smaller beam at this fre- 
quency. The narrow components at -30 and -lOkms non 
observed on our SMA data, may be unrelated interstellar gas, 
or possibly associated with the extended (8 arcmin) emission 
of the planetary nebula. 

' We have applied a conversation factor of 0.8 K/Jy, which takes 
into account the frequency and the size of the beam, to convert the 
SMA line observations to a main beam temperature scale. 
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Fig. 2. '^CO and '^^CO spectra obtained with the JCMT at 
the position RA(2000)= 17''13'"43.9' Dec(2000)=-37°06'1 1". 
Note the presence of three velocity components in each line. 
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Fig. 3. Averaged CO J= 2 ^ 1 spectrum obtained with the 
SMA, integrated over the full spatial extent of the source. The 
dashed lines mark the velocity range of each of the velocity 
components. The stellar velocity, estimated at -31.5 kms"' (cf 
Sec. 6.2), is approximately the boundary between the medium 
and high velocity components. 

Based especially on the averaged spectrum, we define three 
broad velocity components, at ~ -25 km s ~ -40 km s ' 
and a feature/wing at ~ -50 to -60 km s^L These features 
are present in all four spectra, but are best defined at '^CO 
J= 2 — > 1. The same three components are also seen in the 
CO J= 1 -» of Gomez et al. (1989). The low velocity feature 
is much stronger in the SMA integrated spectrum than in the 
J= 2 — > 1 JCMT spectrum, indicating it is located off-centre. It 
is absent from the J= 3 ^ 2 JCMT spectrum. 

4.2. SMA ^2 CO J- 2 ^ 1 and "CO J- 2 ^ 1 maps 

The SMA channel maps allow us to define the precise veloc- 
ity ranges for each of the three broad components. These range 
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Table 2. SMA 1.3 mm dust continuum and CO 2-1 results. Si^t is the integrated flux within the 5cr contour. The FWHM is the 
deconvolved Gaussian size, assuming position angle is zero. 



Tracer 


Coordinates 


^ peak 




FWHM 




(J2000) 


(Jy/beam) 


(Jy) 


(arcsec) 


Cont. 1.3mm 


17:13:44.45 -37:06:11.1 


0.62 


1.26 


6.0 X 4.6 


CO J= 2 -» 1 


17:13:44.47 -37:06:08.0 






11.3 x4.2 




Fig. 4. Superposition of the CO J= 2 — > 1 integrated inten- 
sity map (colour scale) with the '^CO J= 2 — > 1 integrated 
intensity map (contours), both integrated over the full veloc- 
ity range, i.e. from -64 to -15 km s"'. The peak value of the 
'^CO J= 2 ^ 1 emission is 190 K km s ', while for the '^CO 
J= 2 — > 1 emission the contours go from 20 to 90 % in steps of 
10%, with a peak at 78 K km s"^ 



from -64 km s" to -48 km s" for the low velocity compo- 
nent, from -48 km s"' to -30 km s"' for the medium velocity 
component, and from -30 km s ' to -15 km s ' for the high 
velocity component. The relative intensity of these components 
varies from position to position, resulting from a complex ve- 
locity structure in the source. 

The '^CO J= 2 — > 1 transition is optically thinner, and then 
is a better tracer of the mass distribution and kinematics than 
the '^CO J= 2 — > 1 . In the following we then discuss further the 
'^CO J= 2 — > 1 emission. Figure]?] displays the full integrated 
'^CO J= 2 ^ 1 (colour scale) and '^CO J= 2 ^ 1 (contours) 
emission observed with the SMA. We can see that the '^CO 
J= 2 — > 1 emission is more structured and elongated than the 
'^CO J= 2 1 one. Comparing with the continuum emission, 
we see that the line emission is more extended, and peaks three 
arcseconds to the North of the continuum peak (see Table ]2]i. 
The lower contour levels also show extensions to the South and 
East which are not as strong in the continuum. Part of the dif- 
ference may be caused by the free-free contribution to the con- 
tinuum. 



Fig.]5]shows the low velocity interval of the '^CO J= 2 ^ 1 
(a) and '^CO J= 2 — > 1 (b), as well as the medium and high in- 
tervals of the '^CO J= 2 ^ 1 line (b and d). The low velocity 
component (Fig. ]5t and b) differs significantly from one line 
to another. The '^CO shows only one clear component, which 
is offset to the South-East compared to the brighter (and clos- 
est) '^CO component. The missing emission peak in the '^CO 
is likely too weak to have been detected in these observations 
which have a one cr noise level of ^ 2 K km s Its weakness 
indicates a lower opacity for this clump, which is not apparent 
in ^^CO J= 2 — > 1 where both clumps are optically thick. The 
cause for the position shift is less clear, but it may indicate a 
difference between the location of the temperature peak, traced 
by the optically thick '^CO peak, and the column density peak, 
traced by the optically thin '^CO peak. 

The medium velocity component (Fig.]5]:) matches the dark 
lane seen in the Ha map better than the '^CO, and leaves lit- 
tle doubt about its association with the dark lane. We note also 
a slightly curved shape of this component. The high velocity 
component (Fig. ]5}i) is more extended than its '^CO counter- 
part, and is curved in the opposite direction to the medium ve- 
locity component. The lack of evidence for dust extinction as- 
sociated with this strong CO component suggests it is located 
behind the Ha emitting material. The medium and high veloc- 
ity peaks are situated on either side of the continuum peak, but 
shifted North by a few arcseconds. The location and velocity of 
these two components suggests they belong to an expanding, 
toroidal structure, but the northward shift suggests the torus 
may be weaker or incomplete towards the South. 

In the next section we will discuss the mass of this torus 
based on CO observations 

5. Mass determination 

In order to determine the total molecular mass of the torus we 
estimate the excitation temperature and opacity of the CO lines 
using two methods. These both rely on a good knowledge of 
the observed '^CO and '^CO line ratios. 

Since our JCMT '^CO observations only barely resolve the 
torus, even in the J= 3 ^ 2 transition, we assume that the cen- 
tral JCMT spectra contains all the emission from the torus and 
so we can use these spectra alone to determine the total molecu- 
lar mass. However, note that this is not true for the low velocity 
component which is not part of the torus, as this component 
peaks at 7" South-East. 

Table ]3] gives the JCMT main beam integrated intensities 
of the four transitions, ^^CO J= 3 ^ 2, '^CO J= 2 ^ 1, '^CO 
J= 3 — > 2, and '^CO J= 2 — > 1, integrated over the three ve- 
locity components identified in our SMA data. As noted above 
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Fig. 5. Superposition of the Ho- image (colour scale) with CO integrated intensity maps (white contours) obtained with the SMA. 
In each panel the contours go from 20 to 90 % of the emission peak in steps of 10%. a) '^CO J= 2 — > 1 integrated intensity map 
from -64 to ^8 kms"' with a peak value of 47 K kms"'; b) Same as (a) for the '■'CO J= 2 — » 1 transition. The emission peak 
is 26 Kkms c) '^CO J= 2 ^ 1 integrated intensity map from -47 to -30 kms"', the emission peak is 68 Kkms"'; d) '^CO 
J= 2 ^ 1 integrated intensity map from -29 to -15 kms"', the emission peak is 43 Kkms"'. 

Table 3. JCMT '^CO and '^'CO integrated intensities for the three velocity components in T^b scale. The values in parenthesis 
are the one cr standard deviation uncertainties. 



Source Name 


Velcocity range 
(kms-i) 


,'^C0(3-2) 
int 

(Kkms-') 


,'-C0(2-l) 
int 

(Kkms-') 


,'^CO(3-2) 
int 

(Kkms-') 


,'-'CO(2-l) 
int 

(Kkms-') 


low vel. comp. 
medium vel. comp. 
high vel. comp. 


-64 to -48 
-48 to -30 
-30 to -15 


20.0 (2.4) 

56.1 (2.6) 
40.7 (2.4) 


12.5 (1.4) 
27.3 (1.5) 
15.4(1.4) 


1.8(0.3) 
174 (0.3) 
6.2 (0.3) 


2.6 (0.2) 

8.7 (0.2) 
3.0 (0.2) 



the CO J= 3 ^ 2 intensity is likely underestimated for the low 
velocity component. The JCMT beam filling factors can be 
calculated quite accurately for each velocity component us- 
ing the source size measured in our SMA observations. These 
beam filling factors are listed in Table |4] assuming that the 
size of the source is the same for both CO species and both 
transitions. The filling factor corrected line ratios are given in 
Table |5l Also, the comparison of the JCMT and SMA '^CO 
J= 2 — > 1 integrated flux shows that the SMA observations re- 
cover ~ 80% of the JCMT flux, demonstrating once again the 
compact nature of NGC 6302. 



Table 4. Beam filling factors (BFF) 



Source Name 


BFF CO J= 2 ^ 1 


BFF CO J= 3 ^ 2 


low vel. comp. 


0.62 


0.93 


medium vel. comp. 


0.43 


0.64 


high vel. comp. 


0.31 


0.57 



5.1. Analytical estimates of the kinetic temperature 
and opacity 

The ratio of two different transitions of the same molecule 
constrains the excitation temperature, as given by Eq.(3) of 
Levreault (1988), assuming that the beam fiUing factors and 
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Table 5. Integrated line ratios. The values in parenthesis are the one cr standard deviation uncertainties. 



Source Name 


^''C0(3-2) ,^'^C0(2-l) 
int ' inl 


^'-'C0(3-2) ,^'-'C0(2-l) 
inl ' int 


^'^CO(3-2) ,y"CO(3-2) 
inl ' int 


^'-C0(2-l) ,y"C0(2-l) 
int ' int 


low vel. comp. 


1.07 (0.25) 


0.46 (0.11) 


11.11 (3.18) 


4.81 (0.91) 


medium vel. comp. 


1.38 (0.14) 


1.34 (0.05) 


3.22 (0.21) 


3.14(0.24) 


high vel. comp. 


1.35 (0.22) 


1.12 (0.13) 


6.56 (0.70) 


5.13 (0.81) 



the opacities of the transitions are known. For simplicity, we as- 
sume that the excitation temperature is uniform along the line- 
of-sight and that the transitions are thermalized, which allows 
us to infer the opacity of one transition given the opacity of the 
other one. 

For the NGC 6302 case, we know that the '^'CO lines are 
not optically thick (i.e. r < 1). So assuming an opacity of 0.5 
for the '^CO J= 3 — > 2 line we find a kinetic temperature of 

= 25+5 K for the high velocity component and Tk - 35+5 K 
for the medium velocity component. For the low velocity com- 
ponent we find Jk ^ lOK, but because of the missed flux in 
the '^CO J= 3 ^ 2 spectrum (Sec. |5]l we consider this value 
underestimated by a factor of ~ 2 and believe that kinetic tem- 
perature of Tk - 20 K is more realistic for this component. 

With this range of kinetic temperatures inferred from the 
'^CO ratio, it is then possible to infer the opacity of the '^CO 
lines. We find a maximum opacity for the '^CO J= 3 — > 2 
transition of unity for every velocity component, taking into 
account the uncertainties on the observed ratio. 

The opacity and kinetic temperature can also be inferred 
from the line ratio of two isotopic substituted molecules for the 
same transition. This method does not require the knowledge of 
the beam filling factor, if the filling factor is the same for both 
species, but it does require knowledge of the relative abundance 
of the species used. Assuming that the beam filling factor is the 
same for both transitions, that the excitation temperature is the 
same for both transitions, and also that the linewidth is the same 
in both transitions, the ratio of the integrated intensities is given 
by Eq.(3) of Myers et al. (1983). 

Although the terrestrial value of At,, the abundance ra- 
tio X[^^CO]/X['^^CO], is 90, for PNe < 30 (Balser et al. 
2002), and 15 is a lower limit for NGC 6302 (Sec. For 
Ah = 15 the observed '2CO(2-l)/'^CO(2-l) ratios lead to a 
'^CO J= 2 — > 1 line opacity of 3 for the high and low velocity 
components and 5.5 for the medium velocity component, with 
a typical uncertainty of +0.5. The corresponding opacities for 
the '^CO J= 3 -> 2 transitions are then 3 + 0.5 and 7.5 ± 0.5 for 
the high and medium velocity components, respectively. The 
opacity inferred for the low velocity component is ~ 1 + 1.5, 
well below the one inferred from the J= 2 — > 1 transitions. This 
is likely due to the issue with the '^CO J= 3 ^ 2 already men- 
tioned in the previous section. The corresponding opacities in 
the '^CO transitions are thus 15 times smaller, that is ~ 0.1 
for the low velocity component, 0.2 for the high velocity com- 
ponent, and ~ 0.5 for the medium velocity component, with a 
typical uncertainty of +0.05. For an abundance ratio Ai, = 30 
rather than 15, the opacities are approximately doubled. 

The '^CO opacities inferred through this method are 3 to 
15 times higher than the ones inferred using two transistions 



(□) High velocity component (b) Mediurri yelocliy component 




a(2000) q(20Q0) 



Fig. 6. i2cO(2-l)/>3CO(2-l) ratio from the SMA data. The 
star symbol shows the assumed position of the central star 
(Matsuura et al. 2005a). a) High velocity component, b) 
Medium velocity component. 

of the same species. Such differences in the opacity estimates 
lead to a similar dispersion in the mass estimates and indicate 
the need for a more complete model. 

5.2. A two phase model: A hot inner edge and a cold 
outer part 

The most likely explanation for the discrepancy between the 
two analytical methods discussed above is the presence of 
opacity and temperature gradients within the beam (Levreault 
1988). Given the presence of a luminous, hot stellar remnant 
in the centre of NGC 6302, such gradients would not be unex- 
pected. 

Figure |6] shows maps of the ratio '2CO(2-l)/'''CO(2-l) for 
the medium and high velocity components. For both compo- 
nents this ratio is far from uniform over the structure. The high 
velocity component shows higher values of this ratio than does 
the medium velocity component; the peak is located in the 
central part, close to the proposed location of the central star 
(Matsuura et al. 2005a). Direct or indirect energy input by the 
central star appears to be responsible for a temperature gradi- 
ent within the molecular structure, making the gas hotter in the 
inner part and colder in the outer part. 

Figure |6] suggests an inner, marginally optically thin, warm 
component, seen directly through the high velocity component, 
and an outer, optically thicker, cold part, obscuring the warm 
component associated with the medium velocity component. 
This is consistent with the interpretation above that the high ve- 
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locity component is located behind the star, while the medium 
velocity component is in front of star (Sec. 4.2). 

The highest value of the i2cO(2-l)/i3CO(2-l) ratio, found 
for the low velocity component, gives a lower limit to the abun- 
dance ratio of 15. 

To account for the observed non-uniform physical condi- 
tions we have modelled the CO emission assuming it arises 
from two phases with different temperatures. Using the ID ra- 
diative transfer code RADEX (Schoier et al. 2005) we have 
constructed a grid of models varying the kinetic temperature 
and the column density of the CO. The volume density, the 
linewidth and the background temperature were fixed to 1 x 
10^ cm"^, 10 km s"' and 2.73 K, respectively. We used two dif- 
ferent values, 15 and 30, for the abundance ratio between '^CO 
and '^CO. Four grids of models were constructed, one for each 
transition, '^CO J= 2 ^ 1, '^CO J= 3 ^ 2, '^^CO J= 2 ^ 1, 
and '^CO J= 3 — > 2. Each individual model within these grids 
gives the resulting radiation temperature for a medium with 
a different column density and kinetic temperature. The '^CO 
column density was varied from 5 X 10^^ cm"^ to 5 x 10"^ cm"^ 
in steps of 5 x 10'^ cm"^, and the kinetic temperature was var- 
ied from 10 K to 300 K in steps of 5 K. In order to simulate 
a two-phase medium we mixed each point on each grid with 
every other point, adopting a contribution of each phase to the 
observed spectra ranging from 10% to 90% in steps of 10%. 

This procedure generated almost 6x10^ different two com- 
ponent models. We then built for each model, the four line ra- 
tios, i.e. i2c0(3-2)/i2C0(2-l), i3CO(3-2)/i^^CO(2-l), ^^CO(3- 
2)/"CO(3-2), and ^^CO{2- 1 )/i3CO(2- 1). A^^ analysis on these 
grids yields the models which best fit the observed ratios for 
each of the three velocity components. 

Figure [T] shows the results of the radiative transfer mod- 
elling for high velocity component. These plots only shows the 
model ratios fall within the uncertainties of the observed ra- 
tios (Table ID, for the preferred abundance ratio of 15. We can 
clearly see on Fig |7h that two different phases are needed to 
reproduce the observed ratios, one cold phase, with a kinetic 
temperature ~ 20 K and a warmer phase, with a kinetic temper- 
ature ranging from ~80 K to 300 K, depending on the model. 

Figure|7j3 shows the mass distribution of these models. For 
each mass bin we show in red and blue the relative percentage 
of hot and cold material, respectively. Here, hot means a ki- 
netic temperature of Tk > 80 K and cold means Tk < 80 K. The 
most likely mass for the high velocity component is ~ 0.3 
and this mass is largely dominated by hot molecular gas. This 
is consistent with the fact that the high velocity component ex- 
hibits the hot inner part of the torus. 

Figure [8] shows the same plots for the medium velocity 
component. We see that the temperature distribution is quite 
similar to that of the high velocity component, but the mass dis- 
tribution is very different. The most likely mass for the medium 
velocity component is ~ 1.1 M© and this mass is this time 
largely dominated by cold material. This is consistent with our 
picture of NGC 6302 since the medium velocity component 
shows the cold outer part of the torus. Using an abundance of 
30 increases the best fit masses by roughly a factor of 2, but we 
did not find any model for the medium velocity component fit- 



Table 6. Column densities and mass. The mass of the low ve- 
locity component is considered a lower limit. The mass of the 
high velocity components is only likely underestimated 



Component 


(10"^ cm-2) 


(10^' cm-2) 


Mass 
(Mo) 


high velocity 


2.9 


4.4 


0.3 


medium velocity 


10.7 


16.1 


1.1 


low velocity 


0.3 


0.5 


0.06 



ting within the uncertainties of the four line ratios, suggesting 
15 is a better choice value for the abundance ratio. 

The derived mass of the medium velocity component is al- 
most four times higher than the mass of the high velocity com- 
ponent. This can be explained by the optical depth effects: the 
bulk of the cold phase of the high velocity component is hid- 
den by the hot optically thick material in front of it, and thus, 
does not appear in our mass estimate of the high velocity com- 
ponent. However this cold component actually dominates the 
total mass of the toroidal structure. This discrepancy shows the 
limits of our simple two phase model. Although the total mass 
of the molecular content resulting from our models is 1 .4 M©, 
geometrical symmetry suggests that the actual mass is closer 
to 2 Mq, with an estimated uncertainty of 1 Mg given the mass 
distribution of the fit models (Fig. [Tj) and [SJi). This value is 
in very good agreement with the mass estimate from the dust 
continuum emission (see Sect. 3). 

We tried to estimate the mass of the low velocity compo- 
nent in the same way, but excluding the somewhat unreliable 
'^CO J= 3 — > 2 data for this component did not allow us to 
identify a clear mass limit. For this component we have esti- 
mated a lower limit on the mass by assuming a temperature 
Tk = 20 K, LTE, and using the opacity of the '^CO J= 2 ^ 1 
line, ~ 0.35 (obtained in a same way as for the '^CO J= 2 — > 1 
opacity; cf Sect. 5.1). Combining equations (9) and (19) from 
Goldsmith & Langer (1999), and using the '^CO J= 2 ^ 1 in- 
tegrated intensity in Table |3] we find a lower limit of 0.06 Mq 
for the mass of the low velocity component. All these mass es- 
timates are summarized in Table |6] 

6. Kinematics 

The molecular structure surrounding the central star of 
NGC 6302 is elongated in the North-South direction (see Fig.|4] 
andlU and is coincident with the elongated structure previously 
seen in absorption on the Ha image (Matsuura et al. 2005a). 
When we overlay the high and medium velocity components 
observed in CO with the SMA (see Fig.|9t and c), we see that 
these two compoments appear to form the farside and nearside 
parts of a thick disk-like shaped structure surrounding the cen- 
tral star. Based on the angular size observed in '^CO J= 2 ^ 1 
(see Fig.|9j;) we can estimate the physical diameter and thick- 
ness of this structure. Considering that NGC 6302 is located at 
d ^ I kpc (Meaburn et al. 2005), we infer a diameter for the 
disk-like structure of 0.1 1 pc, and a thickness T ^ 0.04 pc. 
However although the morphology is similar to a disk, the ve- 
locity structure is clearly not. 
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High Velocity Component 
Ab=15 



Temperature (,n K) 




High Veiocity Component- 
Ab=15 



Fig. 7. Results of radiative transfer modelling of a two phase 
medium for the high velocity component, a) Temperature dis- 
tribution of the models fitting the four observed line ratios 
within the uncertainties. These models use an abundance ra- 
tio X[i2C0]/X['3C0] of 15. b) Mass disti-ibutions of the mod- 
els fitting the four observed line ratios within the uncertainties. 
The fraction of hot material (i.e. with Tk > 80 K) is represented 
in red, while the fraction of cold material (i.e. with < 80 K) 
is represented in blue. 



6.1. The torus velocity structure 

Considering the medium and high velocity CO components 
as tracing a disk-like structure, the velocity channel marking 
the limit between the two should correspond to the systemic 
velocity of NGC 6302. This velocity is found to be about 
Visi = -30 + 1 kms"'. This is in remarkable agreement with 
the value of -30.4 kms"' derived by Meaburn et al. (2005) 
given the very different methods used. 

For the '^CO J= 2 ^ 1 transition. Fig. |9}3 displays the 
position-velocity (PV) diagram constructed along a North- 
South axis (dashed lines) passing through the center of 
NGC 6302. In this PV diagram we see the presence of a 
ring-shape structure centered at the systemic velocity, Visi- = 
-30 kms"', and at the position 5(2000)=-37"'06'12.5". A sim- 
ilar PV diagram is observed in CO. The low velocity com- 
ponent on the left part of the '^CO J= 2 ^ 1 PV diagrams as a 
fainter peak. The ring-shape structure observed in the PV dia- 
grams does not suggest rotation as it would be the case for a ro- 
tating disk. Compared to the systemic velocity the high velocity 
component, the rear part of the structure, is red shifted, while 
the medium velocity component is blue shifted. This configura- 



Medium Velocity Connponeni 
Ab=15 



(b) 



Medium Velocity Component 
Ab = 15 



in. 



Mass {M,„„) 

Fig. 8. Same as for Fig. |2]but for the medium velocity compo- 
nent. 



tion clearly suggests expanding motion. Rather than a disk, the 
molecular structure seen around the central star of NGC 6302 
is more like an expanding molecular torus. 

For a complete, uniform expanding torus seen edge-on, the 
PV diagram constructed along the torus should ideally reveal 
an ellipse whose extent along the velocity axis directly gives its 
expansion velocity. Actually, the shape and overlap of the front 
and back parts of the torus (see Fig. |9^), as well as the previ- 
ous measurements of the NGC 6302 outflow angle compared 
to the plane of the sky (i.e. 18°; Meaburn et al. 2005) show that 
this torus is seen nearly, but not exactly, edge-on. In Fig.|9l3 we 
plot an ellipse which delimitates the molecular emission. This 
ellipse fits quite well the ring-like shape structure, although we 
note that there is emission missing in the southern part of the 
red shifted emission (i.e. the high velocity component is less 
extended than the medium one; cf Fig.|9^). From these ellipses 
we can infer the mean projected expansion velocity of 9 km s ' . 
Since we know the angle of the outflow axis to the plane of the 
sky, 18°, we can, assuming that the torus axis has the same 
angle, correct the expansion velocity from the projection. But 
since this correction corresponds to only 0.5 kms ' which is 
comparable to the uncertainties, we adopt Vexp = 9 kms~K 
If there was some rotation (or overall velocity gradient) asso- 
ciated with the torus the ellipse in the PV diagram would be 
inclined and possibly twisted. Although difficult to accurately 
constrain without a detailed model. Fig. |9j) indicates that any 
such gradient has a magnitude of < 20 km s"' pc"' . 

The analysis of the lowest velocity CO component is more 
problematic. The emission from this component is less similar 
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Fig. 9. a) Ha map (grey scale) overlaid with the SMA '^CO 
J= 2 1 medium (solid blue lines) and high (dashed red 
lines) components. The black dashed line corresponds to the 
axis along which we constructed the position-velocity diagram 
shown in (b). b) '^CO J= 2 ^ 1 position-velocity diagram con- 
structed along the black dashed axis shown in (a). The ellipse 
delimitates the molecular emission coming from the torus. The 
0" position coiTesponds to (J(2000)=-37"'06'15.9". 



in the '^CO J= 2 -> 1 and '^CO J= 2 ^ 1 transitions than is 
the case for the torus emission. However its greater velocity 
from the systemic velocity of the sources leads us to interpret 
this component as associated with the outflow from the source. 
Interestingly this component is projected onto the source in the 
region where the torus emission is weakest, perhaps suggesting 
this may material being accelerated out of the torus. 

6.2. Position-offset diagrams: ballistic flow 

To uncover the general velocity pattern, we determine for each 
velocity channel the position and flux of the brightest pixel. The 
radial offset of this pixel from the stellar position is measured. 
We assume the star to be located at the position proposed by 
Matsuura et al. (2005a; cf Fig|6|l. This technique reduces dras- 
tically the amount of information contained in the images, but 
has been found to provide a powerful tool, used, e.g., in maser 
studies of circumstellar envelopes. 



We apply this method to the SMA '"^CO J= 2 ^ 1 data, 
which shows best signal to noise ratio. The result is shown 
in Fig. [To] We select all channels with peak flux above 1 
mJy/beam. The left panel shows the positions of these peaks. 
The stellar position is taken as (0,0). Different symbols indicate 
the velocity ranges. Larger symbols are used for stronger peak 
fluxes. The distribution is reasonably well fitted with a tilted 
ring, with a radius of 5.5", an inclination of 18° (Meaburn et 
al. 2005), rotated to a position angle of -5°. 

The right panel of Fig.fTOlshows the velocity versus the an- 
gular separation between the peak and the assumed stellar posi- 
tion. Such diagrams can distinguish between expanding shells 
and polar flows (Zijlstra et al. 2001). Both are present here. The 
solid line is our model for the ring, where the best fit is found 
for a central velocity Visi = -31.5kms"' and an expansion 
velocity Vexp = 8kms These differ slightly from the ear- 
lier fits but are within the uncertainties. The symbols clustering 
around the drawn line shows that part of the CO emission is 
well derscribed by an expanding ring. In the following sections 
we consider the values inferred from the best fit model as the 
reference ones. 

On Fig. [To] (right) we also notice that at higher and lower 
velocity two linear features appear. These show velocity lin- 
early increasing with distance from the star. This is indicative 
of polar flows. The ballistic nature of these can be understood 
as due to the interaction of a slow and a fast wind, where the 
swept-up interface moves at a constant but direction-dependent 
velocity. This gives rise to a ballistic flow pattern (Zijlstra et al. 
2001), seen in bipolar OH/IR stars. 

The velocity range of this component is small in NGC 
6302, compared to typical for bipolar OH/IR stars. The gra- 
dient is ~ 140 km s ' pc But it shows a considerable spatial 
extent, 2-3 times the size of the ring. The linear flows are seen 
on either side of the star, and both approaching and receding 
components with respect to the systemic velocity, are present. 



7. Discussion 

7.1. Comparison with other objects 

Molecular gas is relatively common among planetary nebulae. 
It is seen among the youngest objects where the high density 
material traps the dissociation and ionization fronts within the 
nebulae. As the nebula expands and dilutes, the molecular re- 
gions will be overrun (Huggins et al. 1996, Woods et al. 2005). 
In the most recent and the most sensitive survey, Huggins et 
al. (2005) detect CO emission for 40 nebula out of a sample of 
110, although only a few are as bright as NGC 6302. 

CO torii are known for a number of objects in a similar 
evolutionary phase. M 1-16 (Huggins et al. 2000) has an ex- 
tensive envelope of 0.12 Mq, around a very compact ionized 
region. NGC 7027 also has a high-mass (~ 1 Mq) CO shefl, but 
its ionized nebula is elliptical and lacks the butterfly shape. A 
closer optical morphological analogy is NGC 2346, which is 
also a butterfly nebula: it shows CO distributed in a clumped 
expanding ring with an estimated mass of 0. 1 Mq (Bachiller et 
al. 1989). 
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Fig. 10. Position-offset diagram. Left panel: positions of the peaks. Triangles: -20 > V > -40 km s~^ ; circles: V < -40; (filled) 
squares: V > -20 km s"^ The size of the symbols are a (weak) indicator of peak flux. Only peak flux values > 1 mJy/beam (on 
scale of top panel) are included. The solid line indicates the model for a tilted ring, with a radius of 5.5" an inchnation angle 
of 18° compared to the plane of the sky, and rotated to a position angle of -5°. The (0",0") position corresponds to the location 
where the central star is assumed to be. Right panel: the velocity versus the radial ofi'set to the central star. FiUed symbols are 
positions to the East of the star. The soUd Une is still our best fit torus model where Vlsr = -3 1 .5 km s"' and Vexp = 8 km s"^ . 



During the preceding phase of the post-AGB nebulae, 
where the star has ceased its mass loss and left the AGB, but 
the star is still too cool to cause significant ionization, CO is al- 
ways present. Such objects tend to show strong bipolarity; CO 
velocities up to 200 km s"' are seen (Bujarrabal et al. 2001) al- 
though the spectra are dominated by more sedate material. In 
some cases a Keplerian disk is present (e.g. Bujarrabal et al. 
2005). Woods et al. (2005) show that two distinct classes of 
objects exist in this phase, one rich in difTerent molecules and 
the other only showing a small subset, but both types show CO 
emission. The molecule rich class (objects such as CRL618, 
OH 23 1.8-1-4. 2) shows the highest density torii, similar to that 
in NGC 6302. 

However, NGC 6302 differs from the other objects in the 
class, in that it is more evolved, with a much hotter central star 
and a high mass for the molecular torus. The full line width of 
50 km s"' is typical for planetary nebulae, althought about half 
the line width is due to the bipolar, ballistic flow first detected 
in the SMA data. 



7.2. Torus properties 

The CO emission in NGC 6302 comes predominantly from a 
compact, expanding torus centred on the exciting star. This CO 
emission requires that there are at least two different tempera- 
ture components of the CO. These components are similar in 
temperature to those identified in the dust emission from this 
source by Kemper et al. (2002) suggesting that both dust com- 
ponents are associated with the torus rather than the outflow 
lobes as Kemper et al. suggested for the warmer dust. The prop- 



erties of the CO emission have a number of important implica- 
tions for the origin and evolution of this system. 

The best estimate of the mass of the torus is ~ 2 M© assum- 
ing a '■'CO to '^CO abundance ratio of 15. This mass is similar 
to that inferred from the submillimetre continuum observations 
of Matsuura et al. (2005a), but significantly higher than those 
estimated from previous observations of CO, e.g. Huggins & 
Healy (1989). This discrepancy probably results from a combi- 
nation of previous underestimates of both the integrated inten- 
sity of the CO J= 2 — > 1 line and opacity of the transition. In 
addition to the molecular material traced by the CO, previous 
observations have identified 0.05Mo of ionised HI (Gomez et 
al. 1989), and O.25M0 in low excitation atomic gas traced by 
CII (Castro-Carrizo et al. 2001) (all scaled to the assumed 1 kpc 
distance) associated with the source. The CO torus therefore 
completely dominates the mass of the circumstellar material 
associated with NGC 6302. 

The high mass of the CO torus means that NGC 6302 has 
the highest mass of inner circumstellar material of any PN 
or proto-planetary nebula (PPN) known. The survey of PPNe, 
hkely models for the progenitor to NGC 6302, by Bujarrabal 
et al. (2001) found only two objects with masses of circum- 
stellar comparable or larger than found in NGC 6302, but both 
these objects, AFGL2343 and IRAS-i-10420, are classified as 
yellow hypergiants with broad CO lines without a clear line 
core-wing separation, unlike the PPN sources. However it is 
has been suggested that that the distance to these sources have 
been overestimated, consequently with much lower masses of 
circumstellar material (Josselin & Lebre 2001). 

As mentioned above, the lowest velocity CO component is 
not part of the expanding torus. SpatiaUy it appears close to the 
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region where the torus emission is incomplete, suggestive of it 
being part of the torus which has been distrupted and acceler- 
ated by the outflow from the central star. Evidence that this may 
be occurring comes from the high excitation Unes detected by 
Feibelman (2001). These lines peak at velocities between -54 
km/s and -64 km s ' which would associate them with this 
most blueshifted CO material. Meabum et al. (2005) also finds 
significantly blueshifted forbidden Nn emission at close to the 
position of the lowest velocity CO component. This clump of 
CO material therefore appears to the source of at least some of 
the material seen in the ionised outflow. The velocity of these 
UV and optical fines is quite different from the -23 km s ' to 
-34 km s~^ velocities of the IR recombination Unes and low 
excitation atomic lines (Casassus et al. 2000; Castro-Carrizo et 
al. 2001), which in light of the CO torus would appear kine- 
matically to originate from within the CO torus. 

7.3. Age 

The size and velocity of the CO torus can be used to infer a 
number of timescales associated with the system. The radius 
of the outer edge of the torus is 13", corresponding to 1.3 x 
10^ AU. With the best fit model expansion velocity of 8 km s"' 
this implies that the event which initiated the expansion of the 
torus occurred ~7500 years ago. From the apparent thickness 
of the torus, i.e. 8000 AU, the event which produced the torus 
lasted for about ~4600 years; the event ended ~ 2900 years 
ago. These timescales assume that the limits of the current CO 
emission trace the full extent of the original material of the 
torus and that the apparent inner edge marks the end of the 
ejection event. 

The ballistic flow cannot be dated as its incUnation angle 
with respect to the plane of the sky is not known. If we assume 
it dates from the time after the torus ejection, we find from the 
velocity gradient (140km s"^ pc"^) that i > 30 degrees, where 
i is the angle with the line of sight. 

Meaburn et al. (2005) find a younger age from the expan- 
sion profile of the north-western optical lobe, 1 .7 arcmin from 
the star. They date the lobe to an age of 1900yr. Within the un- 
certainties, this may be similar to the end of the torus ejection 
(as seen in NGC 6537, a similar nebula to NGC 6302: Matsuura 
et al. 2005b). However, the bipolar lobes may also have formed 
later. 

7.4. Progenitor star 

Combining the total mass of the circumstellar material cur- 
rently associated with the nebula, ~ 2 M© , with the estimated 
current mass of the central star, 0.7-0. SM© (Casassus et al. 
2000) gives an absolute lower limit on the mass of the pro- 
genitor of the system of ~ 3Mo. The upper mass limit is 8Mo, 
based on the fact that stars more massive than this do not be- 
come post-AGB stars, but instead ignite their carbon core. This 
mass Umit is consistent with the ~ 4 Mq progenitor required 
to produce the high temperature, high-mass stellar remnant. 
Intriguingly, this means that the bulk of the inital mass of the 
system is stiU present close to the star. 



Bipolar planetary nebulae have on average higher-mass 
progenitors, based on their smaller Galactic scale height of 
bipolar nebulae (Corradi & Schwarz 1995; Stanghellini et al. 
2002). This is consistent with the high derived progenitor mass. 
The formation of bipolar nebulae has been linked to binary 
progenitors (Soker 1998a) and the angular momentum in such 
a system could easily drive the mass loss needed to generate 
the torus. Based on lUE observations, Feibelman (2001) has 
claimed that NGC 6302 is a binary system with a GV sec- 
ondary, however the extensive circumstellar material and the 
extremely high extinction towards the central star makes this 
questionable. The UV flux more Ukely represents scattered 
fight from the central star (see the discussion in Meabum et 
al. 2005). 

An infrared excess seen at the position of the star (Matsuura 
et al. 2005a) possibly represents a very compact circumstellar 
nebula. The excess can be free-free emission (as in Be stars) or 
very hot dust. The presence of circumstellar material located 
within the central cavity can be related to binarity (as in post- 
AGB stars with disks: van Winckel 2002). However, this is in- 
dicative only, and the binary nature of the central star remains 
suspected but unproven. 

Casassus et al. (2000) noted that NGC 6302 does not show 
any evidence for a hot wind or wind blown cavity. This is sur- 
prising given the evidence for the ballistic flow, and the much 
faster flows in the ionized gas (Meaburn et al. 2005). The kine- 
matics of the nebula makes it Ukely that a fast stellar wind did 
exist in the past. 

7.5. Mass Loss and Evolution 

For the progenitor to eject 2 Mq of material over period of 4600 
years (based on the thickness of the torus) implies an aver- 
age mass-loss rate of 5 x 10 '^M© yr concentrated towards 
the equatorial plane. Even adopting the outer radius only of 
the torus to estimate the timescale and reducing the mass to 
mirumum measured, 1.4Mo only reduces the mass-loss rate to 
2 X 10"^ Mo yr~^ This mass loss rate is consistent with those 
measured towards PPNe, although a factor of a few larger than 
typical (Bujarrabal et al. 2001). Interestingly the mass-loss rate 
for NGC 6302 is very similar to that infered from the similarly 
massive, but more spatially extended, CO shell around the hot 
PN NGC 7027 (Masson et al. 1985), alfliough its sheU is much 
less equatorially condensed. 

Models show that mass-loss rates can only approach these 
values during the late AGB stage of evolution. However values 
as high as that seen for NGC 6302 are difficult to produce in the 
models, even for short periods. The NGC 6302 mass-loss rate 
exceeds any value obtained in models by VassUiadis & Wood 
(1993) and is barely reached at the very tip of the distribution 
of the models by Bloecker (1995). 

Mass loss on the AGB is assumed to be driven by radiation 
pressure on the dust. The achievable rates are limited by the 
available momentum in the stellar radiation field: this Umit can 
be exceeded by a factor of two if multiple scattering of photons 
is taking into account (Vassiliadis & Wood 1993). Using the 
lower limit to the mass-loss rate, the expansion velocity derived 
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above and a stellar luminosity L = 10 Lg, we find that the ratio 
between the momentum in the torus and the one available in the 
stellar radiation, is in NGC 6302 : 

^=^=7.8 

Considering that the mass loss is concentrated towards the 
equatorial plane, while the stellar radiation is isotropic, wors- 
ens the discrepancy. We conclude that, despite the low Vexp, 
the mass-loss history of the torus is difiicult to achieve by stan- 
dard radiation-pressure-driven mass loss. A similar argument, 
but with much larger discrepancies, is presented by Bujarrabal 
et al. (2001). On the other hand, the low expansion velocity is 
as expected from a dust-driven wind at very high M (Habing 
et al. 1994). 

In contrast, the energy requirements are relatively easily 
satisfied using a hypothetical binary companion. If we assume 
a 1 Mq companion in a 1 AU orbit, and that the envelope of a 
3 Mq primary star is at the same radius at the onset of mass 
loss, a reduction of the secondary orbital radius to ~ 0.1 AU, 
whilst the mass of the primary reduces from 3 to 1 Mg, will 
suffice to provide the energy to eject the torus. In the final con- 
figuration, the orbital period of the hypothetical binary will be 
P ~ 15 days. We conclude that binary interaction can be a pos- 
sible explanation for the ejection of the torus. A prediction of 
this would be the existence of a companion with an orbital pe- 
riod P ^ I month. 

The only butterfly-type nebula with a known binary nu- 
cleus is the above-mentioned NGC 2346, with a period of 1 6 
days. This object also has a CO ring, albeit with considerable 
lower mass than seen in NGC 6302. It has an unusal period: al- 
most all other known planetary nebulae with close binary cen- 
tral stars show much shorter periods, requiring common en- 
velope evolution (Zijlstra 2006). A period of 16 days requires 
binary interaction (as the final orbit is smaller than the AGB 
star) whilst avoiding common envelope evolution which would 
lead to shorter orbits. However, whether this is common for the 
butterfly nebulae is not known. 

8. Summary 

High angular resolution observations of '^CO J= 2 — > 1 and 
^'CO J= 2 — > 1 have detected a massive torus centred on the 
exciting star. The observations indicate a lower limit of 15 
for the abundance ratio '^CO/'^CO and so by implication, 
'^C/"C. A detailed analysis of the emission in the J= 2 -> 1 
and J= 3 ^ 2 transitions of '^CO and '^CO indicates that the 
emission arises from material at two different temperatures. In 
total the torus contains ~ 2 M© ± 1 Mg of material, a very 
similar mass to that previously derived from observations of 
the dust continuum emission from the source. From our best 
fit model, this torus is expanding with a velocity of 8 km s"^ 
which leads to a dynamical age of the torus of 7500 yr, with 
the iimer edge of the torus being ejected ~ 2900 years ago. This 
ejection may have coincided with the eruptive event proposed 
to have formed the systems bipolar lobes. We speculate that the 
ejection of the inner edge of the torus also terminated the fast 
CO outflow which, while typical in the PPN progenitors of PN, 



is not detected in NGC 6302. The current mass within the PN 
implies that the progenitor of the system had a mass of > 3M0 
and during the formation of the torus sustained a mass loss rate 
of 5 X lO-'^Mgyr-i 

The derived mass-loss history is difficult to fit using mass 
loss driven by radiation pressure. We find that binary interac- 
tion is a plausible mechanism leading to the ejection of the 
torus. 

Ultimately the torus is hkely to become photoionised by the 
central star, as is happening to the CO shell around NGC 7027. 
The photon dominated region where the UV from the star pen- 
etrates the CO torus, which is traced by the warm CO compo- 
nent observed here and the low excitation atomic lines observed 
by Castro-Carrizo et al. (2001), is likely to support a rich and 
complex chemistry, making NGC 6302 an interesting target for 
further molecular line observations. 
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